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1. Introduction 


When a shock wave impinges upon a wall, it penetrates the boun- 
dary layer along the surface and both the shock wave and the boundary 
layer are changed from their undisturbed states. If the boundary layer 
remains unseparated, these mutually induced changes take place in a 
small interaction region. For a turbulent boundary layer, it has been 
established [1-8] that an asymptotic description of the interaction re- 
gion requires a three layer structure. In the outermost layer, com- 
prising most of the boundary layer, pressure forces are much larger 
than forces resulting from Reynolds or viscous stresses so the govern- 
ing equations are those for an inviscid flow. For the limit process 
to be considered, the solutions for this inviscid flow region are those 
given in Part I of this paper [9] , hereafter referred to as (I). Imme- 
diately adjacent to the wall is the wall layer, in which viscous and 
Reynolds stresses dominate to lowest order. Between the? 2 e two layers 
is the Reynolds stress sublayer [l] (referred to as the blending layer 
in reference [2] ) in which momentum transfer toward the wall is carried 
out by turbulent means (Reynolds stresses); the dominant terms in the 
equation of motion are the Reynolds stress, pressure gradient, and 
inertia terms. 

This paper is concerned with the analysis of the flow in the two 
inner layers, the Reynolds stress sublayer and the wall layer, the goal 
being the calculation of the shear stress at the wall in the interaction 
region. As indicated above, the limit processes considered are those used 
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in (I). Thus, if c is equal to the nondimensional difference be- 

tween the velocity and the critical sound speed in the flow external to 
the boundary layer, and u t is the nondimensional friction velocity, we 
consider limit processes such that u t « € « 1. In previous analyses 
for € « u (Reference [ll) and £ = 0(u ) (Reference [2]) it was found 

T T 

that it was not possible to formulate an asymptotic criterion for shock 
induced separation. Here, it will be shown that even for € >> u^ 
there is no apparent asymptotic separation criterion. However, exam- 
ple calculations will be used to show that the equation derived for 
the wall shear stress may be used to predict conditions for incipient 
separation with reasonable accuracy. 

It is worthwhile reiterating the fact pointed out in (I) that for an 
unseparated boundary layer the solutions in the inviscid and inner la- 
yers are uncoupled. Because the inner layers are so thin, the change 
in pressure across them is negligible to the order retained and so the 
solution for the pressure found in the inviscid layer in the limit as 
the wall is approached is indeed the wall pressure. With this pressure 
distribution known, then, solutions in the inner layers may be found, 
leading to a relation for the wall shear stress. Thus, the unseparated 
flow case is a weak interaction problem. This is not the case for a 
laminar boundary layer and occurs for the turbulent boundary layer 
because the wall layer is so thin that the upstream influence of the 
interaction causes negligible lifting of the fluid from the wall; that is, 
to the order retained the V component of the velocity is zero, in the 
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in viscid layer, as the wall is approached. This point will be discussed 
again later. 

In order to complete the formulation of the problem in the inner 
layers, it is necessary to specify a closure condition. Here, we use 
a mixing length model, including . the van Driest damping factor, to 
write an eddy viscosity [10], Such a closure model appears to give sat- 
isfactory results as long as the flow is unseparated [ll] and has the 
virtue of simplicity; when the flow is separated, use of such a model 
gives results which have the correct trends but which do not agree well 
with experiment. 


2. Solutions in the Inner Layers 

As in (I), transonic flow over a flat plate with a turbulent 
boundary layer is considered, with a normal shock wave intersecting 
the boundary layer; an adiabatic wall is assumed as are conditions 
such that the total enthalpy may be taken to be uniform throughout the 
flowfield. Nondimensional Cartesian coordinates X and Y are measured 
parallel and perpendicular to the wa^J. respectively, with the origin at 
the point where the normal shock wave intersects the boundary layer. 
Lengths are made dimensionless with respect to the distance from the 

X 

leading edge of the flat plate to the shock impingement point, L, and 
Cartesian velocity components U and V with respect to the critical 
sound speed in the flow upstream of the shock wave and external to the 


boundary layer (hereafter referred to as the 
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overbars indicate dimensional quantities. The mean temperature, T, 

density, p, pressure, P, and viscosity coefficient, p, are referred 

— ❖ — * 

to their critical values in the external flow, e.g. , T , P , etc. We 

e e 

write the Reynolds number. Re, in the usual fashion and for convenience 

•$* 

define a Reynolds number parameter. Re , as follows 


Re = (^) 

P ' 
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Re = < £ -5r il >. 
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The te 'm <p'V' >/p is included in V, where the primes denote fluc- 
tuating quantities and the bracket denotes the average value. The 
friction velocity, u t , is made dimensionless with respect to a , and 
is defined in terms of the external flow density as follows 
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where t is the shear stress at the wall in the undisturbed flow at 
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X = L, and where c^ is the corresponding skin friction coefficient de- 
fined as shown. Finally, the external flow velocity and Mach number 
are written in terms of a parameter, t , as 


•age ® 
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U = 1 + 6 
e 




(3a) 


1 - (^(U-l) 
V l e 


(3b) 


where for transonic flow, e « 1. As in (I), the problem considered 

here is one for which u_ « e « 1. 

T 

In both inner layers to be considered here, the characteristic 
thickness of the region is small compared to its characteristic length. 

As a result, normal Reynolds and viscous stress terms may be neglec- 
ted compared to the corresponding shear stress terms and the trans- 
verse pressure gradient is negligible, to the order retained in the anal^- 
ysis* The solutions to which these inner layer solutions must match 
in a direction normal to the wall are those solutions found in (I), 
expanded in the limit as y = Y/6 -► 0, where 6 = 6/JL is of the or- 
der of the boundary layer thickness. In the limit as x = X/A -*■ - oo, 
where A = A/L is of the order of the extent of the interaction region, 
the solutions must match with the corresponding relations in the undis- 
turbed boundary layer. It is seen, then, that the flow problem in the 
inner layers of the interaction region is formulated as a boundary la- 
yer problem with a known pressure gradient. This also helps explain 
why an additional layer (Reynolds stress sublayer) is necessary in the 
turbulent boundary layer case. That is, in either the laminar or tur- 
bulent interaction, there is an outer layer in the interaction region 
where pressure forces dominate over shear forces, and inviscid flow 
equations hold to lower orders. Obviously, then^ solutions in the outer 

layer do not satisfy the no-slip condition at the wall and a new boundary 

*>< 

layer must be considered at the wall. In laminar flow, a boundary 
layer is described asymptotically by a single layer and so only one 
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so called viscous sublayer is needed in the laminar interaction. How- 
ever, a turbulent boundary layer has a two layer asymptotic structure; 
as a result, two inner layers are needed to describe this boundary 
layer in the interaction region. The Reynolds stress sublayer is the 
equivalent of the velocity defect layer, as will be seen. 

With the above remarks in mind, it is possible to write a simpli- 
fied set of governing equations in which only those terms needed in 
either of the two layers considered here are retained. They are as 
follows 1 
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(4c) 


T + (*y)U 2 = *f- 


P = pT = P(X) 


Y Re* u ? 

D = { 1 - exp( - — -) } 


(4d) 

(4e) 


(4f) 


where y is the ratio of specific heats, D is the damping factor, and 

2 

kl =0.41 is the von KarmcCn constant. Since terms of order u will 

t 

. ^ 

be retained in the solutions, it should be pointed out that terms such 
as < p 1 U' > /p in the continuity equation (4a), and < p 1 T 1 > in the 
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equation of state, (4e), which are of order u^, are not included because 
perturbat: ns from the undisturbed flow values of each of these terms 
would be ot higher than second order. Since the undisturbed flow solu- 
tions are considered known to second order, and we are interested only 
in the perturbations from the undisturbed flow, i - is not necessary to 
include the terms in question. As mentioned previously, *.t is assumed 
that the wall is adiabatic, and turbulent and laminar Prandtl numbers 
are unity, so that the stagnation enthalpy is constant, as in equation 
(4d). 

As shown in (fl# for the case € / u > > 1 the distance from 
the wall to the sonic line is exponentially small compared to the thick- 
ness of the boundary layer. Sir :e the extent of the upstream influence 
of the interaction region is ordered by the thickness of the subsonic 
region, the upstream influence is confined to a region, hereafter 
referred to as the inner region, which is exponentially small in the x 
direction compared to the main part of the interaction region, hereafter 
referred to as the outer region. That is, in the x direction, the inter- 
action region actually consists of two regions, one thin compared to 
the other; in the y direction, each of these regions is subdivided into 
the three layers mentioned previously. Following the procedure em- 
ployed in (I), the solutions in the Outer region will be shown here 
in some detail. Because the upstream influence is confined to the inner 
region, the flow entering the shock wave in the outer region is simply 
the undisturbed flow at the point ifi question. Inner region solutions. 
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which are found using precisely the same methods employed in the outer re- 
gion, are given in reference [12] . 

Reynolds Stress Sublayer 

In the Reynolds stress sublayer, which is intermediate to the outer 
inviscid flow layer and the wall layer, inertia terms, are balanced by both 
the pressure gradient and Reynolds stress terms in the equation of motion 
in the flow direction. The extent of the outer interaction region is X = 0(A), 
where, as shown in (I), 

A = b 6 6 = 0(u ) (5a, b) 

o t 

b Q = (y+D 1 ^ 1/2 (1 - (2y+l)« /4 + ' ) (5c) 

If the dimensionless (referred to L) thickness of the layer is taken to be 

Y = 0(4), say, then since the Reynolds stress = Ofu^) and from (I), 

d P/3X = 0(u t /A), the fact that the pressure gradient and Reynolds stress 

A 

terms in Eqn. (4b) must be of the same order indicates that 6 = 0 (u t A). 

A A 

Here, for convenience, we define 6, y and x as follows: 

8 = u A ( 6a ) 

T 

Y = 6 y X = Ax (6b, c) 

The solutions to which those in the Reynolds stress sublayer must 
match as y oo are those in the outer inviscid flow layer, written in the 
limit as Y/s - y -*■ 0. The equations for the U component of velocity and 
the pressure, (Eqns. (3.9), (3.12 ), (3.19 ), (3.21) and (4.2) Q f (I) are 
summarized here for convenience. Thus, 
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U(x, 0) = (1 +€ ) 1 + u T [ (1 + 2y 6 + ...)u Q1 (y) 

+ (1 + (y -1) c + . . . ) u^(x) ] + . . . 
- K6 — ln(C fix) + ... 

IT O 
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1 + y (2c + (2y -l)e 2 + . . . ) 

- u^y (1 + (2\ — 1) € + . . . )u^(x) + . . . 

+ K 6 In ( C 5 x) + . . . 

ir o 

f °° u 01 ( ^ )dr 7 
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(7a) 


(7 b) 


(7 c) 


where P =l-\c+0(e)is the dimensionless pressure in the external 
6 

flow. The function u^(y) describes the variation of the velocity from its 
value in the external flow in the velocity defect layer in the undisturbed 
flow. That is, if represents this velocity component, it may be ex- 
panded as 


U a = U e + \ u 01 (y) (8) 

and u^(y) is the variable part of the velocity distribution in the velocity 
defect layer in the corresponding incompressible boundary layer [13] . 

The form given by Coles [14] is used here 


u 0 i(y) = K [In y - II (1 + cos Try)] 0 < y < 1 


(9) 


= 0 


y > i 
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where II is Coles* profile parameter. The lart terms in Eqns. 

(7a) and (7b) are due to the curvature of the wall, i. e., for a wall with 
convex curvature described locally by 

Y = - -“KX 2 + ... (10) 

where K « 1 and K -*■ 0 as u -* 0 and e -► 0 such that K/« -*■ 0. The 

T 

value of the constant C is found from the solution for the flow field exter- 

o 

nal to the boundary layer. 

If Eqn. (7a) is written in terms of the Reynolds stress sublayer varia- 
ble, y, the result is 

A 

U(x,0) = (1 +€)“ 1 +u T [(l + 2y € + ...){~ ln| +ft 01 (^) +...) 

+ (l + (y -i)« + . . . )u^(x)3 + . . . 

- K6 — In (C fix) + ... (11) 

IT O 

where 

O 01 ($) = K _1 (ln y - 211) (12) 

Thus, equation (11) is the equation to which U(x, y) must match as y -*• oo. 

As mentioned previously, 3 P/8 Y = 0 to the order retained here. This is 
easily derived from the equation of motion in the Y direction (e.g. , see 
reference Hvnce, the pressure as written in equation (7b) holds 

throughout both the Reynolds stress sublayer and the thinner wall layer. 

In view of the form of equations (11) and (12), the general expansions 
for U and V are written as follows 
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A 

U(x,y) = (1 +e ) 1 + u t (^ In | + u 01 (y) + u 1 (x,y)] 


+ c U T ^ ln(j)b ^^(x, §) + e y) + • 


+ K6 u lc ( x » y ; 6) + . . . 


(13a) 


V(x, y) = v j v 1 (x, y) + . . . (13b) 

The corresponding expressions for the temperature, T(x, $), and density, 

P (x# y)> are found by substituting equation ( 13a) in the energy equation, (4d) 
and substituting the resulting expression for the temperature and equation 
(7 b) in the equation of state, Eqn. (4e). If the expansions for U, V, P, and 
p and stretched variables x and y are substituted into Equation (4b), the 
governing equations for u^, u^, and u^ are found. Thus, 
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where from Eqn. (7b), P^, and are defined as follows 


Pj(x) = - V u^x) 


P u(x) = - y (2y -lju^x) 


(14b) 


(14c) 


(14d) 
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(lSc) 


P (xj 6) = ^ta(Cfi) 
lc TT O 

A 

It should be noted that both and P. really denote two terms, one of order 

lc lc 

In 6 and one of order 1. They are written as one here for convenience. Also, 

2 * 

it is found [12] that = 0(e u t ), thus confirming the result used in Part I 
that, to the order considered here, V(x, 0) = 0 in the outer inviscid flow 
layer. 


Insofar as &^(x, ■/ ) is concerned, it is seen from Equation (11) that 

-*• Uj(x) as $ -*■ co. It will be shown later that the same functional dependence 

must hold as y -*■ 0. Since P^ = P^(x), the solution which satisfies both 

matching conditions and the governing equation, (14a), is 
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This result has been used in deriving equations (14b) and (14c). It iseasily 
shown [1,12] that the solutions to equations (14b) -(14d) may be written as 
follows 


= 2 W n uT) exp{ '^fe) }de 
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2 x(x-e) 
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where the B^(£ ) are functions to be found by matching. As ^ -*■ oo, the inte- 
gral terms in each of equations (17) go to zero exponentially and it is seen 
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that the remaining terms match with their counterparts in equation (11). 

As x -** 0, for ^ = constant, tho solutions satisfy the shock wave jump con- 
ditions to the relevant order, as they should. As $ -► 0, one finds [l] the 
following asymptotic behavior for the integrals 

x 0 <V 

/ 0 (^iT exp{ -i7(xlT lde ~ - B i (x) + 8 i w + ••• < 18a > 
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where v = Euler* s constant = 0.57721. 

*e 

The solutions for U may thus be found from equation (13a), (16), and 
(17). Since, as mentioned previously, one can find the density and temper- 
ature in terms of the velocity, using the energy equation and the equations 
of state, it is seen that a complete analytical solution may be found for the 
Reynolds stress sublayer in the outer region, valid to terms of order € u t - 
It should be noted that the continuity equation could be used to find the term 
of order in V; since it is not used anywhere in this analysis, the solution 
for V is not included. Finally, it is of interest to write the solution for U 
in the limit as y -*• 0 for later use in matching with the wall layer solution. 
Thus, 


13 


u = (l+o ‘‘ +u T ( K " 1 bi| +6 01 (}) t^w] 

A * 

+ e u t /c _1 ln( | )[2y - B^xjln^) + g^fx) + . • * ] 

+ « u^Iy-Du^x) + 2yft Q1 (J) - B u (x)ln(^-) + g u (x) 

+ K6[- ~ln(C o 6x) - (B ic (x)ln 6 + B lc (x)ln(-^) + g^x) In 6 

g^;x) + (19) 


Wall Layer 

At the wall, Reynolds stresses are zero and the skin friction is, of 
course, due entirely to viscous si ss. Immediately adjacent to the wall, 
then, is a layer in which, as the wall is approached, momentum transfer is 
accomplished less and less by turbulent means and more and more by 
molecular mechanisms* In this layer, Reynolds and viscous stresses 
are of the same order. The flow entering the interaction region in this 
layer has a velocity U = O (u^) and this order holds in the interaction re- 
gion as well. If the thickness of the layer is taken to be Y = 0(6), then by 
equating the orders of the Reynolds and viscous stress terms in equation 
(4b), one can show that 6 = 0[(Re u t ) ]. Here, in order to write 6 in 


terms of familiar quantities, we set 
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With these orders for Y and U, and since in the interaction region X = 0( A ) 
and 9P/9X = 0 (u t /A), it is seen that, even though a pressure gradient 
exists, the only terms in equation (4b) are the Reynolds and viscous stress 
terms, to the order retained. The resulting equation is easily integrated 
to give 


~ ~ a u 2 

pD(y)U y fx) 
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(21b) 


Y = 


6y 
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where, as indicated in Eqn. (21b), the shear stress at the wall, t (x), is 

w 

made dimensionless with its value in the undisturbed boundary layer at 
X = L, so that as x = X/A -* - oo, t -► 1. Equation (21a) is essentially 
the same equation used in references [l] and [2]; the only difference lies 
in the closure conditions used. 

With the orders mentioned above for U, P, X, and Y, and for 
p = 0(1), it is easy to show [12] that V = CKu^fi/A) and to corroborate equa- 
tion (4c) to the order retained. Since U = Ofu^), then from the energy equa- 

2 

tion, (4d), it is seen that T = + Ofu^) and from the equation of state, 

2 

equation (4e), then, that variations in p in the Y direction are also 0 (u t ). 

2 

Hence, to order u^, p = p^ and as pointed out previously by Melnik and 
Grossman [2] and Adamson and Feo [1], the fact that p / p leads to the 
result that limit process expansions in the wall layer do not match with 
corresponding expansions from the Reynolds stress sublayer; thus, this 
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difficulty arises only in compressible flow. The difficulty may be overcome 
by taking advantage of the range of validity of equation (21a). That is, in 
any intermediate limit 6 « Y « 1>, equation (21a) is still the governing 
equation; it is necessary to retain additional terms only for Y = 0$). 

Hence equation (21a) may be used to derive solutions which will match with 
those found using limit process expansions in either limit; Y = 0(6) or 
Y = 0(6). Although the methods of solution used by Adamson and Feo and 
Melnik and Grossman are equivalent, the latter* s method is more straight- 
forward and will be used here. 

It is clear both from physical arguments, and from consideration of 
equation (21a) that as Y increases such that y » 1, the viscous terms be- 
come negligible compared to the Reynolds stress terms; this is borne out 
by using the solution to be derived to compare the two terms. Also, for 
y >> 1, the damping factor D is represented u 'y unity plus exponentially 
small terms (eqn. 5). Finally, the density may oe written in terms of the 
velocity and pressure,in general, by using equations (4e) and (4d). Since 
P = P(x) to the order retained, equation (21a) may be integrated to give. 


L /L? v 1 / 2 /-^ 




U(x, y) = r Sin {p (p- ) 7 (~p / p ) 7 u T (~ In y + B(x))> 


(22a) 


(22b) 


where T^ = (^y 4-1) / 2 is the temperature at the wall and T^/T^ can be 
found in terms of U from the energy equation, (4d). B(x) is a function of 
integration which should be evaluated by matching equation (22a) with the 
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limit process expansion solution valid in the wall layer (y = O0))^ that is, 
it should be found as a result of integration from the wall to the y value in 
question, using the boundary conditions at the wall. However, it is only 
necessary to evaluate B(x) to lowest order here, and this may be done by 
noting that if there were no shock wave, then B(x) = C, the constant from 
the undisturbed flow wall layer solution (eqn. (2.8), in (I)). Here, then, 
B(x) is written in terms of an asymptotic expansion 

B(x) = C + w(u t ) B^(x) + . . . (23) 

where <jj (u ) -*■ 0 as u -*■ 0. 

r T 

Since equations (22a) must match with equation (19), iu the limit as 
y -► oo, $ -*• 0, it is seen that must have an expansion of the following 
form 

A 

2 16 
T (x) = 1 + a. 6 + a_€ + . . . + U T (x) + 6 U ~ ln(c ) T , (x) 

W 1 2 T 1 T K O If 

+ € u t (x) + . . . + K6 t (x) + . . . (24) 

T 11 1c 

That is, as y -*• co, such that u^x^ln y = 0(1), 
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u In y = u In — + u In — + u In y (25) 

t 7 t 6 t ~ t 7 ' 7 

o 

where, as shown in (I) (equation (2.9)) 



(~ ) l/2 r sin^fU /T*} - u (~ + C) 
1 e' t k 


w 


(26) 


and the expansion for shown in equation (24) follows to insure the indica- 
ted matching. Thus, if equation (7b), (23), (24), (25), and (26) are 
substituted into equation (22a) and the resulting equation is compared with 
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equation (19) term by term, the unknown parameters and functions in equa- 
tions (24) and (19) may be found. The resulting solution for t , the shear 
stress at the wall in the outer interaction region, and the corresponding 
values for the B.^ and g^ (calculated once the are known using equation (18b)) 
in equation (19) are 

■ a 1 <V 1) 2 a l 

t (x) = 1 + a. e + ; e + . , . - u ~~r u (x) 

w 12 r 2 1 

A a, 

+ € u t »c Hn (-|)(2v - a x )(y + j - -~) 

si 3ft 

+ € u t | u 1 (x)[~-(l - 3y - ) + y(y + ^)] 

a 

- 2nfc _1 (2y - a A )(v -~j) 

a a 

+ k " 1 (2y+ 1 - ~|)(V - ~ |)(ln x - Y e + In 2k) j + . . . . 

2a x 

+ K 6 In (C fix) + . . . (27a) 

n o 


a l ” - 4 V (sin 1 ~ ) 1 + 2y (27b) 


3 1! 
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(27c) 


a. 


B 11 

= K + 1 - -i ] 

(27d) 

8 11 

= B n (lnx-v e ) 

(27e) 


where u^(x) is given in equation (7c). Equation (27a), then, is the solution for 
T w (x) in the interaction region, including the effects of curvature in the ex- 
ternal flow field. It has, in most respects, the same form as the equation 
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derived by Melnik and Grossman [2]^differing mainly in the order of the various 
terms, the inclusion of specific analytical solutions at each order of approxi- 
mation, and the inclusion of the curvature terms. 

The order to make numerical calculations for a given Reynolds number. 
Re, and external flow Mach number, M^ = 1 + (y+l)e /2 + ...» it is necessary 
to provide relations for u^ and 6 in terms of Re and € . One of the required 
equations is equation (26), with equations (20) for 6; the other is given in (I), 
(eqn. (2. ll), 2. 12)). This equation, with the values of the integrals as given 
by Cebeci and Smith [10] is r jpeated here for completeness. 
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U+nr [i - (u /r) ] 
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(28) 


. c 

where T^/T^ = 1 - (U^/T) from the energy equation, (4d). Finally, it is neces 
sary to write an equation for the viscosity, p(T), to be used in equations (20). 
Here, p = T* 1 was used, with calculations being performed for n =3/4. 

Finally, it should be noted that although the solutions presented here are 
found to orders of approximation such that pressure gradient and inertia terms 
were not retained in the equation of motion in the wall layer, higher order 
solutions involving these terms have been investigated [12]. It was found that 
the first terms to involve the pressure gradient were of order 6/A in U and of 
order 6/u^A in t • Thus, they give very small corrections to the solutions 
presented. 
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3. Numerical Calculations and Separation Criterion 

The variation of t with x for various values of external flow Mach nam- 
w 

ber (and thus c ) and Reynolds nunrb ers representative of mo dern aircraft are 
shown in figure 2. The numerical computations were carried out using equa- 
tions (27a) for p = T^^, y = 1. 4, II = 1/2, C = 5 and K = 0; these values 
seem to be suitable for flow over a flat plate. The curves show the general 
features found experimentally in the interaction region. That is, t^(x) goes 

through a minimum, say (t ) . ; as M increases, (i ) . decreases, 

6 w min e w mm 

while as Re increases (t ) increases. Thus, the effect of increasing M 

w min e 

and therefore the strength of the shock wave is to decrease the value of t 

w 

everywhere in the interaction region and hence to force the flow toward sepa- 
ration; increasing Re gives the opposite effect. It should be noted that 

t -► 1 as x -*• 0 because the solution shown in Eqn. (27a) is for the outer 
w 

interaction region. A solution for valid in the inner region can be written 
in terms of the solutions for the pressure perturbations in the inner in viscid 
flow region [12]; The solution is found in precisely the same manner as 
that illustrated here for the outer region and results in a solution similar to 
that given in Eqn. (27a). Finally, a composite solution for could be writ- 
ten, using the solutions valid in the inner and outer regions. Because of the 
limit processes considered in this work (e » u^), this composite solution 
would show only a small variation in for x < 0. However, since analytical 
solutions cannot be obtained for the pressure in the inner region, no solution 
for in the inner region has been included here; 

It is not possible to compare the solution for with experimental results 
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for a completely two dimensional unseparated flow because none are available. 
In those cases where the flow was apparently unseparated (e.g., references 
[15, 16]), T was not measured, and in more recent work, where t has been 
measured (e. g. , references [17,18,19]) the flow is separated. In separated 
flow, the shock wave takes on a lambda configuration near the boundary layer 
and a relatively strong pressure gradient develops in the Y direction in the flow 
external to the boundary layer [18]; the flow picture is quite different from 
the unseparated flow case considered here. 

Although experimental results for truly two dimensional flow are not 
available for comparison, there is one set of measurements in a tube in which 
the flow is approximately two dimensional [20]. Thus, if R is the dimensionless 
(with respect to L) radius of the tube, 6/R = 0.04 to 0.08j in addition, the 
changes in the core flow (external to the boundary layer) due to the rapid in- 
crease in the boundary layer displacement thickness through the interaction 
region give corrections which are asymptotically of higher order than those 
retained in Eqn. (27a). In presenting the tube data, Gadd fitted power law 
velocity profiles to the measured profiles and inferred values of 6 (dimen- 
sional boundary layer thickness) immediately upstream of the interaction. 

Using equations derived using power law profiles, he also gave Reynolds num- 
bers associated with the tunnel stagnation pressure and Mach number for 
each test. Skin friction measurements were derived from Stanton tube 
measurements. In determining the flow parameters to be used in calculating 
for comparison with each of Gadd' s experiments, it was decided to use 
Gadd' s values of 6, Re, M^, and stagnation pressure, as being a self- 
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consistent set of data to calculate the necessary c , u^, 6, and II for use in 

t . Thus, it is easy to show that if u = T 30 , 
w 


M 


U 


Re 6 = =F 


at 




(1 - 


e. -n-1/2 
2 ' 


a 6 
w 

v 

w 


(29) 


Here, a^ and are the dimensional speed of sound and kinematic viscosity 

respectively, evaluated at the wall temperature (atmospheric temperature, 

taken to be 59°F. ) and in the case of v , at atmospheric sea level density; 

P refers to atmospheric pressure. Using the given values of M , P, 6, and 
at c t 

t , Eqns. (3), (26), (with Eqn. (20) for 6), (28) and (29) were used to calcu- 
late the equivalent t , u^, 6, and II . From equation (2), then, the correspond 
ing c could be calculated. The ratio of the calculated C. to the value in- 
ferred from the standard tube measurements ranged from 1 . 05 to 1. 29 in 
four cases reported by Cadd (Figures 25 to 28, reference [20]). For this 
reason and because of uncertainties in the calibration of the Stanton tube, it 
was decided to compare values of c f /c^, which is equal to r ^ as given in 
equation (27a). The results of this comparison are shown in figure 2, for the 
case = 1.15,, Re - 7 x 10^, (Figure 25, reference [20]). The remaining 
parameters are given under figure (2). The point X'/6 =0, defined b’ Gadd 
as the position at which P w /P = 0. 528, was found by using equation (7b). 

It is seen that the measured upstream influence is not small. That is, 

V^ C fa is not small for X < - 1 say, as required for this theory, so that even 
if the solution for for x < 0 were available, it is not expected that it would 
give good agreement. In fact, using the above mentioned parameters. 


22 


6^/6 = 0.5, where 6^ is the dimensionless distance to the sonic line in the 
undisturbed boundary layer; evidently the values for Re, € , and II do not form 
a good combination for comparison with the theory. On the other hand, a 
slight unsteadiness in the position of the shock wave could have contributed 
to the slow variation of the measured c ^/ upstream of and in the neighbor- 
hood of the minimum. Nevertheless, the value and the position of the minimum 
of c^/ c are predicted quite accurately. Downstream of the minimum the 
comparison is fairly good; in this regard, however, it is interesting to note 
that the negative curvature seen on the calculated curve but not on this partic- 
ular experimental curve, is a feature found in other experimental results 
which could not be used here because small separation bubbles existed. 

It is of interest at this point to consider the problem of predicting condi- 
tions under which the interaction brings the flow to the point of incipient sepa- 
ration. First, it is seen from equation (27a) that there is no asymptotic condi- 
tion for incipient separation; that is, unlike the laminar case, in which 
-1/5 

€ = 0(Re ) is the asymptotic criterion [2l] # there is no relation between e 

and Re which holds in the limit as Re -+■ oo as a condition for separation. 

This is an important difference between the two flows, and if is of interest to 
investigate the reason for its occurrence. The effect of the interaction, 
through the induced adverse pressure gradient, is to slow the fluid. In the 
boundary layer, then, the stream tubes must become wider and, due to the 
coi itraint of the wall, the V velocity component increases at points away 
from the wall, causing the outer flow to lift away from the wall also. In the 
laminar case, the thickening of streamtubes is greatest in the viscous sublayer. 
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The resulting V component of velocity is large enough that the flow external 
to the boundary layer is affected to lowest order so that the external and 
boundary layer flov s must be considered simultaneously, i. e., a strong 
interaction results [21, 22], No matter how large Re becomes, this strong 
interaction occurs, with the thickness of the viscous sublayer and boundary 
layer decreasing as Re increases, according to their asymptotic dependence 
on Re. The sublayer momentum flux and viscous stresses decrease and the 
strength of the shock wave necessary to cause enough displacement of the fluid 
to result in separation decreases as Re increases. In the turbulent case, even 
for € > > u^, the interaction is a weak interaction to lowest order be- 

cause the wall layer is so thin. Thus, until separation occurs, the outward 
displacement of the fluid in the wall layer due to the interaction is too small 
to cause any effect in the lowest order solutions in the flow external to the 
boundary layer. A strong interaction does not occur until a separation bubble 
exists. Since there is no mechanism through which variations in the wall 
layer and external flows can interact, before a separation bubble is 
formed, it appears that no asymptotic criterion exists for incipient separation. 
However, it may be that such a criterion will result from an asymptotic solu- 
tion for the separated flow problem in the limit as the size of the bubble shrinks 
to zero. 

Although the solution for T w ( x ) does not give an asymptotic 

criterion for separation, there remains the possibility that conditions 

for incipient separation can be found simply by assuming that equation (27a) is 

an accurate solution for t (x) at values of M and Re near separation. It is 

w e r 
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clear from figures 1 and 2 that the solution shows the correct form 
with a minimum value, and it is possible to calculate the correspond- 
ing values of Re and M (i. e., e) for which (t ) - 0, the condition 

for incipient separation. It must be emphasized that equation (27a) is not 

being used in an asymptotic sense in such a calculation; thus in order for 

t to go to zero, one or more terms in the expansion must become as large 

as the first term. Instead, we consider equation (27a) as being a good approxi- 

3 2 

mation to t^Ix ) in a numerical sense as long as € and u t (the orders of 

the first terms neglected) are small compared to one. 

To illustrate the use of equation (27a) for t (x) to predict conditions 

v/ 

for separation, we choose the remainder of Gadd' s tubes flow experi- 
ments in which ^ was measured [20], That is, Gadd presented four plots 
of C £ vs. X'/6 (Figures 25-28, reference [20]), the first of which is shown 
in figure 2. In each case he also performed oil -flow experiments, which 
indicated that in one case (used in figure 2) the flow was not separated, but 
that in the three other cases, separation did occur. Although the plotted values 
of p£ did not indicate the occurrence of separation in these three cases, it 
should be noted that the values of c^ were inferred from measurements from 
a Stanton tube aligned facing the flow; thus, accurate reverse flow measure- 
ments could not be made. Calculations of t were made for each of 

w 

these cases, using the same method for calculating the necessary param- 
eters, as mentioned in the discussion of figure 2. The resulting values 

for (t ) = ( c,/c, ) . for each case are as follows! 

w min f fu rmn 
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(1) 

M =1.27 
e 

Re = 10 7 


(t ) . = 

w min 

- 0.080 

(2) 

M =1.26 
e 

Re = 1.2? x 

io 7 

(t ) . = 

w min 

- 0.020 

(3) 

M =1.34 
e 

Re = 1. 93 x 

io 7 

(t ) . = 

w min 

- C 344 


Thus these calculated results indicate that in ail three cases the flow is sepa- 
rated, in agreement with the oil-flow experiments. In case (2), the extent 

of the region where t < 0, i. e. , the extent of the separation bubble, appears 

w 

to be very small; the flow is barely separated. 

If we denote by the Mach number of the external flow at incipient 

separation, equation (27a) may be used, with the condition that ( T w ^ rn ^ n = 

to find M as a function of Re. A typical result is shown in figure (3) for 
es 

K = 0 and II = 1/2, i.e., for conditions associated with flow along a flat 
plate. It is seen that according to this prediction, M increases as Re 
increases. This result is in agreement with measurements made by Roshko 
and Thomke [23] for supersonic flow at high Reynolds numbers. The magni- 
tude of the increase in M over a large range of Re, however, is small 

6 S 

enough that this result could help explain the conclusion that there was little 

or no variation with Re, reached by Settles, Bogdonoff, and Vas [24]. 

The effects of curvature on M , as predicted by equation (27a) can also 

© s 

be compared with experimental results. Evidently, the only data available 
are those presented in figure (37) of reference [20], reproduced here as 
figure 4. The value of the coordinate along Ihe abscissa, t, is given as 
t = 2. 6 6 /R where R is the radius of curvature; in view of equation (10), 
then, one can write t = 2.66K. Although there is no dependence on Reynolds 
number shown, it is assumed here that the range of Reynolds numbers is 
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10 6 to 10 7 and calculated results are given for both values. Finally, on an 
airfoil with supercritical flow, the flow is accelerating up to the shock wave; 
since II depends on the pressure gradient in the undisturbed flow upstream of 
the interaction, the value of II on an airfoil will be different for different 
curvatures. For zero pressure gradient II = 0. 5, whereas for highly accel- 
erating flow II is smaller and can become negative [25]. Therefore, at t = 0, 

II - 0.5, (K = 0) and at t = 0. 015, it was decided to use a value of II for mod- 
erately accelerated flows, H = 0. The values of M at t = C can be found 

cs 

from figure (3). Those at t = 0. 015, for which K = 0.021 at Re = 10^ and 

7 

K = 0.028 at Re = 10 , were calculated, again using equation (27a). The re- 
sults are shown in figure (4). It is seen that at the conditions associated with 
flow over a flat plate (t=0) the calculated M compares very well with the 

6 S 

value given by the line drawn through the experimental data. On the other 
hand, at higher curvature (t = 0.015) the calculated values are considerably 
less than those found experimentally. In reference [20] there was son dis- 
cussion of the fact that criteria for separation might have been too stringent 
in the co' ed surface cases so that, for example, the point through which 
the drawn line passes at M = 1 . 31 perhaps should have been at M =1. 29. 

6 S 6S 

If this were the case and if negative values of TI were called for, the agreement 
at t = 0.015, would be much better. 

The present results for criteria for shock induced incipient separation 
may be compared with theoretical predictions given by Bohning and Zierep [26], 
who postulated a two iayer model for the interaction region and were 
able to calculate an equation for c.^. Two comparisons were made, both for 
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6 6 

flat plates, at Re values of 10 and 5 x 10 ; these Re are in the range of re- 


sults presented in reference [26]. At. Re = 10 , the predicted values of M 
are 1. 24 and 1. 18 and at Re = 5 x 10^ they are 1. 26 and 1. 30, where the 


es 


M calculated by the present method is the first, in each case. Thus, al- 
es 

6 6 

though the two solutions give the same M at some Re between 10 and 5 x 10 , 

6S 

Bohning and Zierep 1 s solution shows a much greater variation of M with Re 

than that shown here. However, the present results appear to be in closer 

agreement with experimental measurements [23, 24] for a related problem. 

The present theory could not be compared with very recent analytical results 

given by Inger [27] , who also used a two layer model, since conditions for 

incipient separation were not presented. 

Although there appears to be no asymptotic criterion for separation in 

the limit as Re -*■ oo (u^ -► 0) and € 0, there remains the possibility that 

there exists a criterion involving a large M as Re -*• oo. Thus, it is neces- 

c 

sary to consider the behavior of for e = 0(1). Based on the present anal- 
ysis, it is seen that for € = 0(1), u^ « 1, the solution for t^(x) would be 
of the following form in ue outer interaction region: 


t (x) = t (c ) + u t (x;e ) + . . . 
w w , x w, 

d 1 


(30) 


where t (e ) is the value which t (x) approaches far downstream of the shock 
d 

wave. Since the lowest order solutions for the velocities would be of the same 


form in each of the layers as for € « 1, 


it is seen that, from matching 


solutions in the limit u 0, one would obtain 
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W 


r sin {(■ 


V P e 


1/2 -1 

I ' ■ n 


u 


-in ( -jr )} = U d 


(31) 
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Here and are the values of P and U immediately downstream of the 

shock wave in the inviscid flow external to the boundary layer respectively, 

and are thus the values which P and U approach as x -* oo in the interaction 

region. If the jump conditions across a shock wave are used to write and 

P d /P e in terms of and these expressions are substituted into equation (31), 

one obtains an expression for t in terms of U . Thus, 
c v/ , e 

d o 


w 


- 


-1,1 " 

1 +(^)(U 2 - 1) 
c e 


sin (pjj-) 
e 

i - (^)(u 2 - 1) 


sin (~) 

2 e 



(32) 


If U = 1 + € is substituted into equation (32) and the resulting equation ex- 

2 

panded for c « 1, it is found that = 1 + a^ + a^fa^-l)* /2 + in 

d 

agreement with the first three terms of equation (27a). 

It is clear from equation (32) that, since 1 < U < F for 1 < M < co 

© € 

(see eqn. (36)), t f 0 for any M ; instead t goes through a minimum 
W d 6 W d 

value of 0. 512 at M^ = 2. 55, for y = 1.4, and then begins to rise with in- 
creasing M . Hence, there is apparently no limiting value for M as 
6 © S 

Re -»• oo. Moreover, since U -*■ V asM -» oo, it is clear that t -*• ao. 

e e w, 

d 

Recalling the definition of t (x), one can see that this limit means that for 

w 

high Mach number flow the shear stress far downstream of the interaction 
must be large compared to that of the undisturbed flow. This apparent 
anomaly can be explained by considering equation (21a). The density, in the 
first term, can be written as P/T through the equation of state. Now, in the 
wall layer, the temperature differs from the constant temperature of the 
wall by only higher order terms. On the other hand, since 3 P/9 Y = 0 
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through the wall layer and Reynolds stress sublayer, P is the pressure from 

the inviscid flow layer and so varies from P e to P^. As -* oo, P^ /P g oo 

and so from equation (21a), r -*00 also. In general, since the Reynolds 

W d 

shear stress, - p < U* V' >, includes the density, this result appears to be 

independent of the specific closure condition as long as < U' V* > does not 

go to zero as M -*• 00 , and is another significant departure from, the laminar 

case. Experimental verification of the large values of t at high Mach num- 

W d 

bers is given in measurements by Marvin, et al. [28] . 


4. Concluding Remarks 

The use of asymptotic methods of analysis results in a relatively simple 

relation for the shear stress at the wall in the interaction region. This relation 

may be used to predict conditions for incipient separation. In order to obtain 

the proper variation of t vs. x, which includes a minimum in t , it is neces- 

sary to include terms of higher order than the first approximation; evidently 

this would be the case also if one were to calculate r for the case u = O U ) 

w T 

[ 2 ] • 

Although the range of parameters in available experiments does not allow 
for exhaustive testing of the theory, comparisons which could be made are en- 
couraging; more accurate results should be obtained at the high Reynolds 
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numbers associated with modern transonic '-'r w raft. 
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Figure Titles 


Figure 1 


Figure 2 


Figure 3 


t vs x for various values of Me and Re, for flow over a flat 

W 

plate. (Eqn. (27a) with II = 0,5, K = 0). 

Comparison of calculated (eqn. (27a)) and experimental (ref. 
[203) values of vs X'/f>. Experimental conditions; 

Me = 1.15, Re = 7 x 10^, wall temperature = 15°C (59° F), 

P te = .r 137 « 9kP a g au 8 e ( 20 P si g)» 6 = 0.305 cm (0.12 in). 

Corre sponding calculated parameters; * ■= 0. 120, u t = 0.03968, 
6 = 0.01634, n = 0.312, K = 0. X»/6 = 0 at P w /P te = 0.528 
from eqn. (7a). 

M es (Mach number at incipient separation) vs Re for flow 
over a flat plate, calculated using eqn. (27a) with ( T w ) rtl £ n =0, 
n = 0.5, K = 0. 


Figure 4 Effect of curvature parameter, t, on Mach number for incipient 
separation, M eg , from Reference [20] . Triangles and circles 
correspond to different airfoils. Present calculations shown 
as follows: at t = 0, ■ - Re = 10^, n = 0. 5, K = 0; x - Re - 10^, 
n = 0. 5, K = 0. At t = 0.015, •- Re = I0 6 , n = 0, K = 0.21; 
x - Re = 10 7 , n = 0, K = 0.28. 
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Figure 1. t w vs x for various values of Me and Re, for flow over a flat plate. 
(Eqn. (27a) with IT = 0.5, K = 0). 





